NACA TN 2642 


loan only, case F I L 

COPY 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 2642 


APPLICATION OF LINEAR ANALYSIS TO AN EXPERIMENTAL 

INVESTIGATION OF A TURBOJET ENGINE WITH 

PROPORTIONAL SPEED CONTROL 

By Marcel Dandois and David Novik 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 

Febriiary 1952 


2319 


lA 


NATIONAL ADVTSORT COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2642 

APPLICATION OP LINEAR ANALYSIS TO AN EXPERIMENTAL 
INVESTIGATION OF A TURBOJET ENGINE WITH 
PROPORTIONAL SPEED CONTROL 
By Marcel Dandols and David NotTc 


SUMMARY 

Results of an analytical and a sea-level experimental investigation 
of a turtojet engine with, proportional speed control are presented. 

Linear analysis and description of the engine as a first-order linear 
system, proved adequate for analytical prediction of the response and the 
stahility of the controlled engine, althou^ instahllity calculations 
were found to he much more critical than response calculations. 

On the basis of a compromise between speed of response and oscilla- 
tions, an optimum loop gain was found. Increased loop gain increased the 
speed of response and decreased the speed error but ultimately led to 
instability characterized by an essentially constant frequency and con- 
stant amplitude oscillation. Operation near the limits of stability 
required a decrease in control gain with decreasing engine speed. 


INTRODUCTION 

'A current approach to synthesis of control systems for turbojet 
engines Involves the techniques of linear analysis and synthesis 
described in references 1 and 2 and their bibliographies. This 
approach is predicated upon a knowledge of the dynamic characteristics 
of the various components of the controlled system. Accordingly, con- 
siderable effort has been expended in obtaining the dynamic characteris- 
tics of gas-turbine engines from experimental and analytical studies 
(references 3 to 6) . These investigations of gas-turbine engine dynamics 
indicate that the engines may be considered linear for substantial excur- 
sions from steady-state equilibrium points and that within the limits of 
experimental accuracy the engines appear to be first-order systems. 

The adequacy of the linear techniques and data discussed in the 
references may be ascertained by comparison of analytical results based 
on these techniques and actual experimental data. An experimental 
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investigation of a turtojet engine controlled "by various types of control 
system has been initiated at the NACA Lewis laboratory to provide experi- 
mental data for this comparison and simultaneously to determine the char- 
acteristics of the different types of control system utilized. 

The results obtained from the initial phase of the experimental pro- 
gram, a sea-level experimental Investigation of proportional speed con- 
trol, are presented herein. Proportional speed control was selected as 
the initial phase of the project bebause engine speed is the controlled 
variable in most engine control systems now in use and because propor- 
tional control is a dominant factor in virtually every type of control 
system. The results obtained from this investigation should be applicable, 
at least in part, to studies of any specific control system for turbojet 
engines . 

The use of a purely proportional control consisting only of an ampli- 
fication factor would not contribute any fundamental change to the dyna- 
mics of the engine and if the engine is considered as a first-order system 
instability could not be obtained and examined. In order to introduce the 
possibility of instability, the control system was permitted to contribute 
several lags to the controlled engine system in addition to the purely 
proportional control element. A nonlinearity other than the engine was 
also included so that the adequacy of linear analysis was severely tested. 


APPARATUS 

Engine 

The turbojet engine used in this investigation was of the axial -flow 
type having a twelve-stage compressor, a s ingle -annuliis combustion cham- 
ber, and a two-stage turbine. Experimentally determined steady -state 
performance curves of the engine are shown in figure 1. 


Control System 

The engine control system consisted essentially of an electronic 
amplifier with variable gain for amplification of an error signal and an 
oil system with a pressure regulating servovalve positioned by the out- 
put voltage of the electronic amplifier, such that a given pressure 
exists for each output voltage. The control oil pressure acts on the 
servopiston of a motor-driven constant -pressure fuel pump (reference 7) 
which supplies a proportional fuel pressure to the engine. For conven- 
ience of analysis and data presentation, the oil system is considered as 
part of the fuel system and therefore the electronic amplifier can be con- 
sidered as the proportional control element. A schematic drawing of the 
controlled engine system is shown in figure 2. Lags in the system Include 
the oil system, the fuel system, and the ELC (resistance-inductance- 
capacitance) filter, as well as the engine. 
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Ins t rument at ion 

Be c order . - All transient data were recorded on a multichannel 
recording oscillograph. The oscillograph elements had a sensitivity of 
1-lnch deflection per 0.43 millivolt and were damped to an effective 
time constant of approximately 0.005 second. 

Engine speed . - Transients in engine speed were measured hy record- 
ing the output of a small tachometer generator on the recording oscill- 
ograph. The high frequency ripple from the generator was filtered out 
of the speed measuring circuit with an EC (resistance-capacitance) cir- 
cuit filter. In the closed loop (fig. 2 ), where engine speed is fed 
hack into the loop and where phase shift and amplitude attenuation become 
Important, an ELC filter was utilized. The frequency response character- 
istics of both filters are shown in figure 3. 

Pressures . - Transients in fuel pressure and control oil pressure 
were measured with bridge-type strain-gage pressure pickups energized by 
a 3000 cycle carrier system. The apparent time constant of the pressure 
measuring circuits was 0.02 second. 


ENGIKE DYNAMICS 

Engine dynamics were obtained .by observing the response of the engine 
variables to approximate step inputs in fuel pressure and to sinusoidal 
fuel pressure changes of variable frequency. These teats were run with 
an open loop system, that is, engine speed feedback was not used and the 
electronic control was omitted. 

Typical data showing the response of the uncontrolled engine are 
shown in figure 4. Variables other than control oil pressure, fuel 
pressiire, and engine speed were obtained for further analysis but are 
not discussed in this report. From these data, the corrected time con- 
stant of engine speed with respect to fuel pressure was obtained and is 
shown in figure 5 as a function of final corrected engine speed. The 
time constants shown were obtained from the frequency response charac- 
teristics of the variables as indicated by harmonic analyses of the data 
obtained with digital computer machines and by the response to sinusoi- 
dal inputs. Frequency response characteristics of the engine, corrected 
for the filter characteristics, are shown in figure 6 for three engine 
speeds. Frequency response characteristics of the fuel system, including 
the oil syston, were obtained from sinusoidal inputs and are shown in 
figure 7 for three engine speeds. 
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RESULTS Al® DISCUSSION 
Effect of Loop Gain on Response and Error 

Response. - The controlled engine system shovm schematically in fig- 
ure 2 was operated over a range of proportional control settings (gain 
settings) for several conditions of engine speed. Transients were induced 
hy means of a step input in voltage (Ng of fig. 2) into the electronic 
control and the responses of the variahles were recorded. Data typical 
of the responses obtained are shown in figure 8. Traces shown for speed 
error Ng-N also represent the negative of engine speed Inasmuch as the 
speed setting Ng is constant after the initial step transient. 

The loop gains shown in the data are defined as the product of the 
gains at zero frequency of all the components in the closed loop (fig. 2) 
such that : 


(control (oil (fuel 

amplifier) system) pump) 


( engine) (tachometer ) (tachometer 

amplifier) 


loop gain 


4olts\ /ib/aq. in.\ 4b /sq ln.\ ( rpm ) 
Uolt j I volt j Ub/sq in.y (ib/sq in^ 

dimensionless q.uantity 


^Yolts^ 
i rpm j 


/volts 1 
I VO It ; 


These loop gains are calculated from step input Ng data from the 
expression 

^ initial error 
loop gain = ‘ 1 


( 1 ) 


This equation is valid for all controlled systems that do not contain an 
Integral term because the final error in such systems is zero. 

In order to show the comparative effect of different loop gains on 
the speed of response of the controlled engine, an arbitrary indication 
of the Initial speed of response was selected rather than the time 
required for attainment of equilibrium. This selection was made because 
the time required for equilibrium to be reached does not necessarily 
decrease with increasing loop gain, Inasmuch as the oscillations result- 
ing from an underdamped system at higi gain may take longer to damp out 
as the gain is increased. The arbitrary criterion of response time was 
taken as the time required for the engine speed to first attain 63 per- 
cent of its ultimate equilibrium speed and may be compared with the time 
constant of the uncontrolled engine. The effect of Increasing loop gain 
on the response time of the controlled engine (time required for engine 
to attain 63 perceiit of equilibrium speed) is shown in figure 9 for three 
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engine speeds. Figure 9 indicates that the controlled engine responds 
more rapidly as the loop gain increases, althou^ the improTement in 
response hecomes less pronounced at the hi^er values of loop gain. The 
response of the controlled engine at 84 percent maximum engine speed is 
faster thau at 90 percent maximum engine speed "because the time constant 
of the engine increased sli^tly as the speed was increased above 84 per- 
cent maximum engine speed. 


The effect of lags and hi^er-order systems on response is shown hy 
a comparison of the experimental data at 84 percent maximum engine speed 
with the theoretical response of a first-oid.er system (fig. 10). The 
first -03Tder curve is "based only on the engine time constant at 84 per- 
cent maximum engine speed and all other components are neglected. This 
curve was calculated from the operational equation 


N / ^ 


1 +TiP 


1 + 


1 + T^iP 


( 2 ) 


where is the engine time constant at 84 percent of maximum engine 

speed (approximately 1.2 sec). 

The comparison of the experimental data with the theoretical first- 
order curve indicates that the response time Increases for hl^er-order 
systems and that the improvement in response with Increased gain hecomes 
more marginal for a system with lags . 

The response of a theoretical second -order system, the curve of which 
is shown in figure 10, is calculated from the following operational 
equation ; 

G 

(1 + Tip)(l + Tgp) 

(1 + t^3^p)( 1 + “TgP) 



^ (P) 


’where ^ again ^ 

engine t±me constant 

Tg first-order approximation of oil and fuel systems obtained from 
63 -percent point of response of fuel pressure to a step input 
(0.28 sec, fig. 4(h)) 
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The hi^er-order system of the actual controlled engine has been 
simplified analytically by considering it as only a second-order system. 
The excellent agreement between the experimental data obtained and the 
calculated second-order curye Indicates that many simplifying assumptions 
may be made in calculating response without causing appreciable deviations 
from the experimentally obtained response of a complex physical system. 


Becatise the system becomes more and more underdamped as the loop 
gain increases (fig. 8), the magnitude and duration of the oscillations 
at hi^ gain tend to negate the advantages of the improved Initial 
response. A criterion is therefore req_uired that will compromise between 
the fast initial response at hi^ gain and the oscillations incurred. 

Such a criterion may reasonably be based on the time integral of the 
square of the speed error produced by a given disturbance. The premise 
is made that the best compromise between response and oscillations will 
occur when this Integral is a minimum. Mathematically^ this relation 
may be expressed as 






dt = minimum 


The speed error is squared in order to give a positive value to the 
Integral when the oscillations result in negative values of engine speed 
with reference to the final engine speed (reference 8) . 


Values of this Integral plotted against loop gain for three engine 
speeds and a small pictorial representation of the meaning of the integral 
are shown in figure 11. For the specific system used in this investiga- 
tion, the values of the integral tend toward a minimum which may be taken 
to occur at a loop gain of about 8. According to figure 9, the region of 
this minimum also corresponds to loop gains at which the improvement in 
response with increasing gain may be considered marginal. 


Figures 9 and 11 indicate that a proportional control setting that 
would give a loop gain of about 8 would result in optimum proportional 
control on the basis of a compromise between initial response and mini- 
mum oscillation. For systems with less lag, this optimum gain would be 
slightly hi^er. 

Error. - An obvious characteristic of proportional speed control is 
that as the gain is increased, the engine speed will more closely approach 
the speed setting, inasmuch as a smaller error will be required 

to support the engine speed at a given value . This characteristic is 
relatively unimportant when the desirability of arriving at an engine 
speed close to the speed setting is considered, because the speed setting 
throttle can always be calibrated or adjusted with allowance for the speed 
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error. However, If the effects of disturbances ether than those initia- 
ted from the control lever are considered, such as those from a plugged 
fuel nozzle or large change in airplane angle of attack, the relation 
between gain and error assumes greater importance. Equilibrium condi- 
tions more closely corresponding to the original conditions will be 
restored and the response will be faster if the gain is high because 
the entire loop is mor.fe sensitive at high gain and only a relatively 
small change in error results in a comparatively large corrective signal 
to the engine. 

Data typical of the engine response to a disturbance other than a 
change in control setting are shown in figure 12, These data were 
obtained by closing a solenoid valve in one of two fuel lines in parallel. 
Because the fuel pump was a constant pressure pump, the increased resis- 
tance to flow resulted in a decreased fuel flow from the pump. 

The data presented show the sudden decrease in fuel pressure followed 
by its gradual increase as the control system attempted to restore the 
original equilibrium conditions . The engine speed following the fuel flow 
disturbance responds more rapidly and more nearly returns to the original 
value for the setting of higher gain. A plot of the difference between 
initial and final engine speeds for vario-us values of loop gain is shown 
in figure 13 for a constant value of disturbance in the fuel system equi- 
valent to a 5 percent maximum engine speed change in the uncontrolled 
engine. The calculated curve was obtained through the use of equation (l). 

Figures 12 and 13 show that as the gain is increased the change in 
error due to random distiurbances decreases . This change in error is the 
difference between the initial and final speeds as plotted on figure 13 
and is actually a change in engine speed only, inasmuch as the speed set- 
ting Ng remains fixed. 


Instability 

Occurrence and natirre of instability . - As the gain settings were 
increased, instability was eventually attained at each engine speed. 
Typical instability data are shown in figure 14. 

Controlled engine instability manifested itself as an essentially 
constant frequency and constant amplitude oscillation. The sinusoidal 
characteristic of this instability is Indicative of the system's linear- 
ity, whereas the fact that the amplitude remains approximately constant 
is Indicative of limits and other nonlinearities within the system. As 
an ideal linear system would oscillate with ever increasing amplitude, 
the presence of nonlinearities has apparently stabilized the system to a 
condition of steady oscillations. 
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Tile oscillations illustrated in figure 14 were obtained at approxi- 
mately the lowest possible loop gains at which instability could be 
attained. At successively hi^er engine speeds, the frequency of oscil- 
lations increased and the magnitude decreased from about 2 to 0.5 percent 
of maximum engine speed. Further increase in loop gain beyond the gain 
at which instability first occurred resulted in a sli^tly increased 
amplitude and frequency of oscillation (these data are not shown) . All 
oscillations obtained represented unfavorable conditions of engine 
operation. 

Effect of increased phase shift in system . - According to linear 
theory, instability occurs when a loop gain of at least unity exists at 
the frequency corresponding to a total phase shift of 180° between input 
(error) and output (speed) of a closed loop system (references 1 and 2). 

In order to indicate the effect of phase lag on instability, the total 
phase shift throu^ the system therefore was increased by substituting 
the large phase lag EG filter for the ELC filter normally used (fig. 3). 
Instability data shown in figure 15 were obtained. 

As a result of the increase in phase lag, the maximum allowable 
loop gain prior to attainment of instability was considerably reduced 
at all engine speeds. A comparison of the loop gains at instability 
for the original system (with ELC filter) and the system with increased 
phase lag (EC filter) is shown in figure 16 over a range of engine speeds. 

Comparison of figure 14 with figure 15 shows that at the same engine 
speeds, the original system oscillated at a hi^er frequency than the 
system with Increased phase lag. The greater allowable loop gain prior 
to instability and the hi^er frequency of oscillation agrees with linear 
instability theory in that reduction of phase shift may be expected to 
displace the 180° phase lag point to correspond to a higher frequency 
and lower amplitude at a hi^er loop gain. 

Prediction of loop gain at instability . - Althou^ the type of 
instability obtained in this investigation exhibited certain nonlinear 
characteristics (see control oil trace, fig. 15(a)), the linearity 
apparent from the data indicated that reasonable estimates of the loop 
gains corresponding to instability could be made. The method of calcu- 
lation consisted in cascading the transfer functions of the components 
within the closed loop system, as represented by their frequency response 
characteristics. By multiplying the amplitude ratios and adding the 
phase lags of the fuel system, the engine, and the ELC filter, the over- 
all frequency response of the system was obtained for each specific 
engine speed. 

The amplitude ratios of each component and of the closed loop system 
were taken as unity at zero frequency. The calculated gain at instabil- 
ity was therefore the factor required to increase the amplitude ratio of 
the closed loop system at 180° phase lag to a value of unity. 
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The freq[uency response characteristics of the components are pre- 
sented In figures 3 (ELC filter), 6, and 7. The cascaded frequency 
response and calculated gain at Instability for the Instabilities of 
figure 14 are shown In figure 17. A comparison of the calculated and 
experimental data appears In the following table: 


Engine speed 
(percent 
maximum) 

Measured 
loop gain at 
Instability 

Calculated 
loop gain at 
Instability 

Measured 
frequency at 
Instability 
( radians / b ec ) 

Calculated 
freq^uency at 
instability 

64 

12.2 

10.0 

4.3 

3.7 

80 

10.0 

12.0 

7.5 

6.4 

88 

9.6 

12.5 

9.8 

7.5 


Comparison of the calculated and actual gains at Instability shows 
that despite nonllnearltles within the physical system, the loop gain 
at Instability can be reasonably predicted from a knowledge of component 
dynamics. Frequency response data for all the components In the system 
were utilized for the calculations, however, and the polar plots of fig- 
ure 17 show that small discrepancies In phase angles could result In 
large changes In the calculated loop gain at Instability. Therefore the 
linear analysis used In the Instability calculations Is adequate If suf- 
ficient frequency response data for every con5)onent In the system are 
available. This conclusion constrasts with the fact that response time 
can be calculated quite accurately despite rather broad simplifying 
assumptions as discussed previously. An example of the difference 
between data required for Instability and response calculations Is that 
Instability would not have been obtained at any value of loop gain for 
the second-order approximation used successfully In the calculation of 
response. 


Maximum Allowable Control Gain 

The variation of loop gain with engine speed limited by Instability 
Is shown In figure 16. For both systems shown the loop gain at Instabil- 
ity decreased with Increasing engine speed. This trend, however. Is not 
directly Indicative of the required variation In control gain setting 
(which may be considered as the Independent variable with respect to loop 
gain adjustment) because of the large variation of engine gain and fuel 
system gain with engine speed (fig. 18). 

If lines of constant control gain are added to figure 16, the shape 
of these curves will reflect the change In engine and fuel system gain 
with engine speed as shown In flgizre 19. At a constant control gain, 
deceleration from a high speed to a low speed could result In a trans- 
ition from stable to unstable operation. For example. If operation Is 
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initially asaimed at 92 percent mazimum engine speed and a control gain 
of 0.2, a deceleration at constant control gain would result in instabil- 
ity at about 80 percent maximum engine speed because the engine gain and 
therefore the loop gain increases as the speed decreases. In order to 
operate near maximum loop gain at all engine speeds, the control gain 
must be adjusted so as to decrease with decreasing engine speeds, a 
trend opposite to that exhibited by the loop gain. 

Wxmbers shown on figure 19 for values of control gain are not quan- 
titatively significant because they Include a lumping of the tachometer 
generator gain, tachometer generator output amplifier gain, and electronic 
control gain. 


Mathematical Representation of Engine 

The correspondence of the frequency response characteristics of the 
engine with theoretical first-order systems may be seen from the phase 
against amplitude plots of figure 6. Deviations are apparent in the 
middle region of phase shift but at phase shifts above 80° the discrep- 
ancies are negligible. Inasmuch as instability occurred at engine phase 
shifts greater than 80°, Instability calculations based on theoretical 
first-order curves would have given the same results at those obtained 
from actual data. First-order linear representation of the engine there- 
fore appears adequate for estimation of controlled engine characteristics, 
including prediction of the loop gain necessary to produce engine 
instability. 


SUMMARY OF RESULTS 

Results from an analytical and experimental investigation of a 
tirrbojet engine with proportional speed control were presented and dis- 
cussed. Description of the engine as a first-order linear system and 
linear analysis based primarily on frequency response techniques proved 
adequate for calculation of response and instability despite certain 
nonlinearities in the controlled engine system. Whereas broad simpli- 
fying assumptions were permissible in calculations of response, the pre- 
dictions of instability required accurate knowledge of the frequency 
response characteristics of every component in the system. 

An optimum loop gain, based on a compromise between speed of response 
and oscillations, was found to exist, and for the specific system investi- 
gated, occurred at a loop gain of about 8. At hl^er loop gains, instabil- 
ity characterized by an essentially constant frequency and constant ampli- 
tude oscillation was attained. The loop gains at which instability 
occurred decreased with increasing engine speed. Because of large changes 
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In engine gain it wae, nevertheless, necessary to decrease the control 
gain with decreasing engine speed in order to permit stable operation. 

As expected, the effect of increasing gain resulted in faster 
response, more overshoot and oscillations, and more satisfactory restor- 
ation of eq.uilihfivmi for a random disturbance. Increased phase lag in 
the system decreased the permissible loop gain limited by instability 
and reduced the improvement of increased gain on response. 


Lewis Fli^t Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 14, 1951 
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APPEKDIX - SYMBOLS 

The follorwlng symbols are used in this report; 

C capacitance 

G loop gain 

L inductance 

N engine speed 

Ng engine speed setting 

N-t engine speed at time t 

final engine speed 

E resistance 

t time 

t^ final time 

t^ initial time 

T temperature 

air flow 

cl 

Wj. fuel flov 

5 pressure correction, ratio of ambient pressure to sea-level 

pressure 

9 temperature correction, ratio of ambient temperature to sea- 

level temperature 

T time constant 
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32 48 64 80 96 

Corrected engine speed, percent maximum 

(a) Variation of corrected fuel flow with 
corrected engine speed. 

Figure 1. - Steady- state performance curves of test engine. 
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(c) Variation of corrected turbine outlet gas 
temperature vith corrected engine speed. 

Figure 1. - Concluded. Steady-state performance 
curves of test engine. 
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Figure 2. - Schematic diagram of closed loop controlled engine systems 
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(t) Phase shift against frequency. 
Figixre 3. - Frequency response of filters. 
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(a) Engine speed, 64 to 68 percent maximum. 



Cb) Engine speed, 80 to 84 percent maximum. 



Control oil 
sure 


Fuel pressure 






4 5 

Time, sec 



(c) Engine speed, 88 to 92 percent maximum. 

Figure 4. - Response of uncontrolled engine to approximate step in fuel flow. 
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Figure 5. - Variation of engine time constant with engine speed. 
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Figure 6. - Frequency response of engine speed to fuel pressure. 
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Frequency, radlans/sec 



-100° -90° -80° -70° -60^ 

Phase angle, deg 


-50° 


(c) Engine speed, 90 percent maximum. 

Figure 6. - Concluded. Frequency response of engine speed to fuel pressure. 
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Figure 7. - Frequency response of fuel pressure to speed setting. 
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(b) Engine speed, 84 percent maximum. 

Flgvu?e 7. - Continued. Frequency response of fuel pressure to speed setting. 
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Figure 8. 



(a) Loop gain, 2. 



Time, sec 
(c) Loop gain, 10. 


- Response of controlled engine to step change in speed setting. Final engine 
speed, 84 percent maximum. 
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Loop gain, G 

Figure 9. - Effect of loop gain on response time. 



Figure 10. - Comparison of calculated response with experimental 
response for various loop gains. 



Integral of speed error squared, 
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Figure 11. - Integral of speed error squared as response criterion for 

various loop gains. 
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If 

Control oil pressure 




("b) Loop gain, 0.51. 



Time, sec 
(c) Loop gain, 2.4. 

Figure 12. - Response of controlled engine to fuel flow disturbance. 
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Loop gain, G 


Figure 13. - Variation of speed error with loop gain following 
step disturbance in fuel system equivalent to 5 percent 
maximum engine speed. 
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Fuel pr»*r.uiire 
Control oil pressure 


Speed err* r 


Time, sec 

(c) Engine speed, 88 percent maximum. 
Figure 14. - Controlled engine instability. 
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(aj Engine speed, 64 percent maximum. 



Figure 15. - Controlled engine instability . System revised by increased phase shift. 
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Figure 16. - Loop gain at instability for various engine speeds. 
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•Figure 17. - Prequen' y response of all loop components In ca^scade 


NACA TN 2642 


35 




to 

cvj 





Phase shift, deg Amplitude ratio 


36 


NACA TN 2642 





Figure 17. - Concluded. Frequency response of all loop components In cascade. 
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Figure 18. - Variation of combined engine and fuel 
system gain with engine speed. 
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Figure 19. - Variation of loop gain with engine speed for constant 

Loop gain 


values of control gain = 


Engine gain x Fuel system gain 
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